extern int dlabad_(doublereal *small, doublereal *large);
